Identical sequence patterns in the ends of exons and introns of human protein-coding genes

Intron splicing is one of the most important steps involved in the maturation process of a pre-mRNA. Although the sequence profiles around the splice sites have been studied extensively, the levels of sequence identity between the exonic sequences preceding the donor sites and the intronic sequences preceding the acceptor sites has not been examined as thoroughly. In this study we investigated identity patterns between the last 15 nucleotides of the exonic sequence preceding the 5' splice site and the intronic sequence preceding the 3' splice site in a set of human protein-coding genes that do not exhibit intron retention. We found that almost 60% of consecutive exons and introns in human protein-coding genes share at least two identical nucleotides at their 3' ends and, on average, the sequence identity length is 2.47 nucleotides. Based on our findings we conclude that the 3' ends of exons and introns tend to have longer identical sequences within a gene than when being taken from different genes. Our results hold even if the pairs are non-consecutive in the transcription order. (C) 2012 Elsevier Ltd. All rights reserved.

Identification of unannotated exons of low abundance transcripts in Drosophila melanogaster and cloning of a new serine protease gene upregulated upon injury

Abstract Background The sequencing of the D.melanogaster genome revealed an unexpected small number of genes (~ 14,000) indicating that mechanisms acting on generation of transcript diversity must have played a major role in the evolution of complex metazoans. Among the most extensively used mechanisms that accounts for this diversity is alternative splicing. It is estimated that over 40% of Drosophila protein-coding genes contain one or more alternative exons. A recent transcription map of the Drosophila embryogenesis indicates that 30% of the transcribed regions are unannotated, and that 1/3 of this is estimated as missed or alternative exons of previously characterized protein-coding genes. Therefore, the identification of the variety of expressed transcripts depends on experimental data for its final validation and is continuously being performed using different approaches. We applied the Open Reading Frame Expressed Sequence Tags (ORESTES) methodology, which is capable of generating cDNA data from the central portion of rare transcripts, in order to investigate the presence of hitherto unnanotated regions of Drosophila transcriptome. Results Bioinformatic analysis of 1,303 Drosophila ORESTES clusters identified 68 sequences derived from unannotated regions in the current Drosophila genome version (4.3). Of these, a set of 38 was analysed by polyA+ northern blot hybridization, validating 17 (50%) new exons of low abundance transcripts. For one of these ESTs, we obtained the cDNA encompassing the complete coding sequence of a new serine protease, named SP212. The SP212 gene is part of a serine protease gene cluster located in the chromosome region 88A12-B1. This cluster includes the predicted genes CG9631, CG9649 and CG31326, which were previously identified as up-regulated after immune challenges in genomic-scale microarray analysis. In agreement with the proposal that this locus is co-regulated in response to microorganisms infection, we show here that SP212 is also up-regulated upon injury. Conclusion Using the ORESTES methodology we identified 17 novel exons from low abundance Drosophila transcripts, and through a PCR approach the complete CDS of one of these transcripts was defined. Our results show that the computational identification and manual inspection are not sufficient to annotate a genome in the absence of experimentally derived data.

Assessing Pathogenicity for Novel Mutation/Sequence Variants: The Value of Healthy Older Individuals

Improvement in DNA technology is increasingly revealing unexpected/unknown mutations in healthy persons and generating anxiety due to their still unknown health consequences. We report a 44-year-old healthy father of a 10-year-old daughter with bilateral coloboma and hearing loss, but without muscle weakness, in whom a whole-genome CGH revealed a deletion of exons 38-44 in the dystrophin gene. This mutation was inherited from her asymptomatic father, who was further clinically and molecularly evaluated for prognosis and genetic counseling (GC). This deletion was never identified by us in 982 Duchenne/Becker patients. To assess whether the present case represents a rare case of non-penetrance, and aiming to obtain more information for prognosis and GC, we suggested that healthy older relatives submit their DNA for analysis, to which several complied. Mutation analysis revealed that his mother, brother, and 56-year-old maternal uncle also carry the 38-44 deletion, suggesting it an unlikely cause of muscle weakness. Genome sequencing will disclose mutations and variants whose health impact are still unknown, raising important problems in interpreting results, defining prognosis, and discussing GC. We suggest that, in addition to family history, keeping the DNA of older relatives could be very informative, in particular for those interested in having their genome sequenced.

PIK3CA exon 20 mutations are associated with poor prognosis in breast cancer patients

OBJECTIVES: The phosphatidylinositol 3-kinase/AKT axis is an important cell-signaling pathway that mediates cell proliferation and survival, two biological processes that regulate malignant cell growth. The phosphatidylinositol 3-kinase CA gene encodes the p110 alpha subunit of the phosphatidylinositol 3-kinase protein. There are phosphatidylinositol 3-kinase CA mutations in several types of human tumors, and they are frequently observed in breast cancer. However, these mutations have not been investigated in Brazilian breast cancer patients. METHODS: PCR-SSCP and direct DNA sequencing were performed to identify phosphatidylinositol 3-kinaseCA exon 9 and exon 20 mutations in 86 patients with sporadic breast cancer. The relationships between PIK3CA mutations and patient clinicopathological characteristics and survival were analyzed. The presence of the TP53 mutation was also examined. RESULTS: Twenty-three (27%) of the 86 primary breast tumors contained PIK3CA mutations. In exons 9 and 20, we identified the hotspot mutations E542K, E545K, and H1047R, and we identified two new missense mutations (I1022V and L1028S) and one nonsense (R992X) mutation. Phosphatidylinositol 3-kinase CA exon 20 mutations were associated with poor overall survival and TP53 gene mutations. CONCLUSIONS: Phosphatidylinositol 3-kinase CA mutations are common in tumors in Brazilian breast cancer patients, and phosphatidylinositol 3-kinase CA and TP53 mutations are not mutually exclusive. Phosphatidylinositol 3-kinase CA exon 20 mutations are associated with poor survival, and they may be useful biomarkers for identifying breast cancer patients with aggressive tumors and for predicting the response to treatment with PI3K pathway inhibitors.

Aberrant transcript produced by a splice donor site deletion in the TECTA gene is associated with autosomal dominant deafness in a Brazilian family

We ascertained a Brazilian family with nine individuals affected by autosomal dominant nonsyndromic sensorineural hearing loss. The bilateral hearing loss affected mainly mid-high frequencies, was apparently stable with an early onset. Microsatellites close to the DFNA8/DFNA12 locus, which harbors the TECTA gene, showed significant multipoint lod scores (32) close to marker D11S4107. Sequencing of the exons and exon-intron boundaries of the TECTA gene in one affected subject revealed the deletion c.5383 + 5delGTGA in the 5' end of intron 16, that includes the last two bases of the donor splice site consensus sequence. This mutation segregates with deafness within the family. To date, 33 different TECTA mutations associated with autossomal dominant hearing loss have been described. Among them is the mutation reported herein, first described by Hildebrand et al. (2011) in a UK family. The audioprofiles from the UK and Brazilian families were similar. In order to investigate the transcripts produced by the mutated allele, we performed cDNA analysis of a lymphoblastoid cell line from an affected heterozygote with the c.5383 + 5delGTGA and a noncarrier from the same family. The analysis allowed us to identify an aberrant transcript with skipping of exon 16, without affecting the reading frame. One of the dominant TECTA mutations already described, a synonymous substitution in exon 16 (c.5331 G<A), was also shown to affect splicing resulting in an aberrant transcript lacking exon 16. Despite the difference in the DNA level, both the synonymous substitution in exon 16 (c.5331 G<A) and the mutation described herein affect splicing of exon 16, leading to its skipping. At the protein level they would have the same effect, an in-frame deletion of 37 amino-acids (p.S1758Y/G1759_N1795del) probably leading to an impaired function of the ZP domain. Thus, like the TECTA missense mutations associated with dominant hearing loss, the c5383 + 5delGTGA mutation does not have an inactivating effect on the protein. (C) 2012 Elsevier B.V. All rights reserved.

High-Performance Liquid Chromatography Under Partially Denaturing Conditions (dHPLC) is a Fast and Cost-Effective Method for Screening Molecular Defects: Four Novel Mutations Found in X-Linked Chronic Granulomatous Disease

Implementing precise techniques in routine diagnosis of chronic granulomatous disease (CGD), which expedite the screening of molecular defects, may be critical for a quick assumption of patient prognosis. This study compared the efficacy of single-strand conformation polymorphism analysis (SSCP) and high-performance liquid chromatography under partially denaturing conditions (dHPLC) for screening mutations in CGD patients. We selected 10 male CGD patients with a clinical history of severe recurrent infections and abnormal respiratory burst function. gDNA, mRNA and cDNA samples were prepared by standard methods. CYBB exons were amplified by PCR and screened by SSCP or dHPLC. Abnormal DNA fragments were sequenced to reveal the nature of the mutations. The SSCP and dHPLC methods showed DNA abnormalities, respectively, in 55% and 100% of the cases. Sequencing of the abnormal DNA samples confirmed mutations in all cases. Four novel mutations in CYBB were identified which were picked up only by the dHPLC screening (c.904 insC, c.141+5 g>t, c.553 T>C, and c.665 A>T). This work highlights the relevance of dHPLC, a sensitive, fast, reliable and cost-effective method for screening mutations in CGD, which in combination with functional assays assessing the phagocyte respiratory burst will contribute to expedite the definitive diagnosis of X-linked CGD, direct treatment, genetic counselling and to have a clear assumption of the prognosis. This strategy is especially suitable for developing countries.

Li-Fraumeni-like syndrome associated with a large BRCA1 intragenic deletion

Background: Li-Fraumeni (LFS) and Li-Fraumeni-like (LFL) syndromes are associated to germline TP53 mutations, and are characterized by the development of central nervous system tumors, sarcomas, adrenocortical carcinomas, and other early-onset tumors. Due to the high frequency of breast cancer in LFS/LFL families, these syndromes clinically overlap with hereditary breast cancer (HBC). Germline point mutations in BRCA1, BRCA2, and TP53 genes are associated with high risk of breast cancer. Large rearrangements involving these genes are also implicated in the HBC phenotype. Methods: We have screened DNA copy number changes by MLPA on BRCA1, BRCA2, and TP53 genes in 23 breast cancer patients with a clinical diagnosis consistent with LFS/LFL; most of these families also met the clinical criteria for other HBC syndromes. Results: We found no DNA copy number alterations in the BRCA2 and TP53 genes, but we detected in one patient a 36.4 Kb BRCA1 microdeletion, confirmed and further mapped by array-CGH, encompassing exons 9-19. Breakpoints sequencing analysis suggests that this rearrangement was mediated by flanking Alu sequences. Conclusion: This is the first description of a germline intragenic BRCA1 deletion in a breast cancer patient with a family history consistent with both LFL and HBC syndromes. Our results show that large rearrangements in these known cancer predisposition genes occur, but are not a frequent cause of cancer susceptibility.

Using the Hard, Randy, and Rittler Test to Evaluate Color Vision in Capuchins (Cebus libidinosus)

The identification of color vision types in primates is fundamental to understanding the evolution and biological function of color perception. The Hard, Randy, and Rittler (HRR) pseudoisochromatic test categorizes human color vision types successfully. Here we provide an experimental setup to employ HRR in a nonhuman primate, the capuchin (Cebus libidinosus), a platyrrhine with polymorphic color vision. The HRR test consists of plates with a matrix composed of gray circles that vary in size and brightness. Differently colored circles form a geometric shape (X, O, or Delta) that is discriminated visually from the gray background pattern. The ability to identify these shapes determines the type of dyschromatopsy (deficiency in color vision). We tested six capuchins in their own cages under natural sunlight. The subjects chose between two HRR plates in each trial: one with the gray pattern only and the other with a colored shape, presented on the left or right side at random. We presented the test 40 times and calculated the 95 % confidence limits for chance performance based on the binomial test. We also genotyped all subjects for exons 3 and 5 of the X-linked opsin genes. The HRR test diagnosed two subjects as protan dichromats (missing or defective L-cone), three as deutan dichromats (missing or defective M-cone), and one female as trichromat. Genetic analysis supported the behavioral data for all subjects. These findings show that the HRR test can be applied to diagnose color vision in nonhuman primates.

Non-classical HLA-E gene variability in Brazilians: a nearly invariable locus surrounded by the most variable genes in the human genome

The non-classical human leukocyte antigen (HLA) class I genes present a very low rate of variation. So far, only 10 HLA-E alleles encoding three proteins have been described, but only two are frequently found in worldwide populations. Because of its historical background, Brazilians are very suitable for population genetic studies. Therefore, 104 bone marrow donors from Brazil were evaluated for HLA-E exons 14. Seven variation sites were found, including two known single nucleotide polymorphisms (SNPs) at positions +424 and +756 and five new SNPs at positions +170 (intron 1), +1294 (intron 3), +1625, +1645 and +1857 (exon 4). Haplotyping analysis did show eight haplotypes, three of them known as E*01:01:01, E*01:03:01 and E*01:03:02:01 and five HLA-E new alleles that carry the new variation sites. The HLA-E*01:01:01 allele was the predominant haplotype (62.50%), followed by E*01:03:02:01 (24.52%). Selective neutrality tests have disclosed an interesting pattern of selective pressures in which balancing selection is probably shaping allele frequency distributions at an SNP at exon 3 (codon 107), sequence diversity at exon 4 and the non-coding regions is facing significant purifying pressure. Even in an admixed population such as the Brazilian one, the HLA-E locus is very conserved, presenting few polymorphic SNPs in the coding region.

Hyperinsulinism/hyperammonemia (HI/HA) syndrome due to a mutation in the glutamate dehydrogenase gene

The hyperinsulinism/hyperammonemia (HI/HA) syndrome is a rare autosomal dominant disease manifested by hypoglycemic symptoms triggered by fasting or high-protein meals, and by elevated serum ammonia. HI/HA is the second most common cause of hyperinsulinemic hypoglycemia of infancy, and it is caused by activating mutations in GLUD1, the gene that encodes mitochondrial enzyme glutamate dehydrogenase (GDH). Biochemical evaluation, as well as direct sequencing of exons and exon-intron boundary regions of the GLUD1 gene, were performed in a 6-year old female patient presenting fasting hypoglycemia and hyperammonemia. The patient was found to be heterozygous for one de novo missense mutation (c.1491A>G; p.Il497Met) previously reported in a Japanese patient. Treatment with diazoxide 100 mg/day promoted complete resolution of the hypoglycemic episodes. Arq Bras Endocrinol Metab. 2012;56(8):485-9

Footprints of a trypanosomatid RNA world: pre-small subunit rRNA processing by spliced leader addition trans-splicing

The addition of a capped mini-exon [spliced leader (SL)] through trans-splicing is essential for the maturation of RNA polymerase (pol) II-transcribed polycistronic pre-mRNAs in all members of the Trypanosomatidae family. This process is an inter-molecular splicing reaction that follows the same basic rules of cis-splicing reactions. In this study, we demonstrated that mini-exons were added to precursor ribosomal RNA (pre-rRNA) are transcribed by RNA pol I, including the 5' external transcribed spacer (ETS) region. Additionally, we detected the SL-5' ETS molecule using three distinct methods and located the acceptor site between two known 5' ETS rRNA processing sites (A' and A1) in four different trypanosomatids. Moreover, we detected a polyadenylated 5' ETS upstream of the trans-splicing acceptor site, which also occurs in pre-mRNA trans-splicing. After treatment with an indirect trans-splicing inhibitor (sinefungin), we observed SL-5' ETS decay. However, treatment with 5-fluorouracil (a precursor of RNA synthesis that inhibits the degradation of pre-rRNA) led to the accumulation of SL-5' ETS, suggesting that the molecule may play a role in rRNA degradation. The detection of trans-splicing in these molecules may indicate broad RNA-joining properties, regardless of the polymerase used for transcription.

Evolutionary history of exon shuffling

Exon shuffling has been characterized as one of the major evolutionary forces shaping both the genome and the proteome of eukaryotes. This mechanism was particularly important in the creation of multidomain proteins during animal evolution, bringing a number of functional genetic novelties. Here, genome information from a variety of eukaryotic species was used to address several issues related to the evolutionary history of exon shuffling. By comparing all protein sequences within each species, we were able to characterize exon shuffling signatures throughout metazoans. Intron phase (the position of the intron regarding the codon) and exon symmetry (the pattern of flanking introns for a given exon or block of adjacent exons) were features used to evaluate exon shuffling. We confirmed previous observations that exon shuffling mediated by phase 1 introns (1-1 exon shuffling) is the predominant kind in multicellular animals. Evidence is provided that such pattern was achieved since the early steps of animal evolution, supported by a detectable presence of 1-1 shuffling units in Trichoplax adhaerens and a considerable prevalence of them in Nematostella vectensis. In contrast, Monosiga brevicollis, one of the closest relatives of metazoans, and Arabidopsis thaliana, showed no evidence of 1-1 exon or domain shuffling above what it would be expected by chance. Instead, exon shuffling events are less abundant and predominantly mediated by phase 0 introns (0-0 exon shuffling) in those non-metazoan species. Moreover, an intermediate pattern of 1-1 and 0-0 exon shuffling was observed for the placozoan T. adhaerens, a primitive animal. Finally, characterization of flanking intron phases around domain borders allowed us to identify a common set of symmetric 1-1 domains that have been shuffled throughout the metazoan lineage.

A new portal for the analysis of human splicing variants

Understanding alternative splicing is crucial to elucidate the mechanisms behind several biological phenomena, including diseases. The huge amount of expressed sequences available nowadays represents an opportunity and a challenge to catalog and display alternative splicing events (ASEs). Although several groups have faced this challenge with relative success, we still lack a computational tool that uses a simple and straightforward method to retrieve, name and present ASEs. Here we present SPLOOCE, a portal for the analysis of human splicing variants. SPLOOCE uses a method based on regular expressions for retrieval of ASEs. We propose a simple syntax that is able to capture the complexity of ASEs.

Distribution of QPY and RAH haplotypes of granzyme B gene in distinct Brazilian populations

Introduction: The cytolysis mediated by granules is one of the most important effector functions of cytotoxic T lymphocytes and natural killer cells. Recently, three single nucleotide polymorphisms (SNPs) were identified at exons 2, 3, and 5 of the granzyme B gene, resulting in a haplotype in which three amino acids of mature protein Q48P88Y245 are changed to R48A88H245, which leads to loss of cytotoxic activity of the protein. In this study, we evaluated the frequency of these polymorphisms in Brazilian populations. Methods: We evaluated the frequency of these polymorphisms in Brazilian ethnic groups (white, Afro-Brazilian, and Asian) by sequencing these regions. Results: The allelic and genotypic frequencies of SNP 2364A/G at exon 2 in Afro-Brazilian individuals (42.3% and 17.3%) were significantly higher when compared with those in whites and Asians (p < 0.0001 and p = 0.0007, respectively). The polymorphisms 2933C/G and 4243C/T also were more frequent in Afro-Brazilians but without any significant difference regarding the other groups. The Afro-Brazilian group presented greater diversity of haplotypes, and the RAH haplotype seemed to be more frequent in this group (25%), followed by the whites (20.7%) and by the Asians (11.9%), similar to the frequency presented in the literature. Conclusions: There is a higher frequency of polymorphisms in Afro-Brazilians, and the RAH haplotype was more frequent in these individuals. We believe that further studies should aim to investigate the correlation of this haplotype with diseases related to immunity mediated by cytotoxic lymphocytes, and if this correlation is confirmed, novel treatment strategies might be elaborated.

New mutations in the GLA gene in Brazilian families with Fabry disease

Fabry disease (FD) is an X-linked inborn error of glycosphingolipid catabolism that results from mutations in the alpha-galactosidase A (GLA) gene. Evaluating the enzymatic activity in male individuals usually performs the diagnosis of the disease, but in female carriers the diagnosis based only on enzyme assays is often inconclusive. In this work, we analyzed 568 individuals from 102 families with suspect of FD. Overall, 51 families presented 38 alterations in the GLA gene, among which 19 were not previously reported in literature. The alterations included 17 missense mutations, 7 nonsense mutations, 7 deletions, 6 insertions and 1 in the splice site. Six alterations (R112C, R118C, R220X, R227X, R342Q and R356W) occurred at CpG dinucleotides. Five mutations not previously described in the literature (A156D, K237X, A292V, I317S, c.1177_1178insG) were correlated with low GLA enzyme activity and with prediction of molecular damages. From the 13 deletions and insertions, 7 occurred in exons 6 or 7 (54%) and 11 led to the formation of a stop codon. The present study highlights the detection of new genomic alterations in the GLA gene in the Brazilian population, facilitating the selection of patients for recombinant enzyme-replacement trials and offering the possibility to perform prenatal diagnosis. Journal of Human Genetics (2012) 57, 347-351; doi:10.1038/jhg.2012.32; published online 3 May 2012